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The complete set of matrix elements for all polarization configurations in top quark 
decays is presented including higher order QCD corrections. The analysis is done 
in the framework of the helicity formalism. The results can be used in a variety of 
circumstances, e.g. in the experimental analysis of top quark production and decay 
at Tevatron, LHC and NLC. Relations to LEP1 and LEP2 physics arc pointed out. 



1 Introduction 

Since its discovery in 1995 the top quark has been an object of increasing 
interest. The production process for top quarks has been analyzed in various 
theoretical studies both for proton and e + e~ collisions. Early references on 
the lowest order cross section areEl for pp collsions andcl for e + e~ annihilation. 
Higher order corrections to the cross section (total cross section, px distribution 
etc.) have been calculated by several groups, Ho for pp andErQ for e + e~. These 
total cross sections do not involve information on the top quark polarization. 
Such spin effects only come in, if one studies distributions of top quark decay 
products. In some cases, spin effect ji. have been studicd_i.e. the distribution 
of the top quark spin vector, in ref. 13 for pp and in ref. □ for e + e~ collisions. 
These latter studies have not been extended to higher orders yet. However, an 
interesting step in this direction has been taken in ref. B, where a Monte Carlo 
progam including final state spin terms has been written. Unfortunately, that 
paper is in fact concerned with higher order corrections to t + t~ production 
in e + e~ annihilation at lower energies and does not take into account axial 
vector couplings. More general results, including axial vector couplings, can 
be found inla, but not in the form of a Monte Carlo program. 

Nothing is known about higher order spin corrections to the processes 
which induce top quark production in proton collisions (light quark annihila- 
tion qq — ► ti and gluon-gluon fusion qq — > tt). In contrast to e + e _ annihilation, 
there are many higher order diagrams involved in these processes, so that the 
calculations are very difficult. 
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Of course, in studying distributions of top quark decay products one has 
to take into account the higher order top quark decay matrix elements, too. 
The spin averaged matrix elements have been calculated including higher order 
corrections. For example, oneloop QCD corrections are known to decrease the 
total width of the top quark by about 8%. If one studies distributions of top 
quark decay products with reference to the production process, one needs in 
addition the spin correlations in the top quark decay matrix elements. In lowest 
order, these spin correlations are known, but in higher order only rudimentary 
results exist EJO. In those calculations, higher order corrections to special 
decay product distributions were calculated. This corresponds to certain linear 
combinations of the polarized matrix elements. The present article describes 
how all the polarized matrix elements can be obtained in a complete and 
systematic way in the framework of the so-called hclicity formalism. 

Within the Standard Model, all couplings of the top quark to other par- 
ticles are completely fixed by its mass and by a few quantum numbers. For 
example, the coupling of the top quark to gluons is a pure vector coupling with 
strength g s , the coupling to the W-boson is a V-A coupling etc. The coupling 
to the W-boson is particularly interesting for this article, because it induces 
the decay t — ► W + b. b-quark mass terms 0(mb/m t ) can very probably be 
neglected in higher order corrections to the decay process, because they are 
known to be small (~ 1%) in leading order. In this approximation it can be 
shown that the b-quark is always left-handed in the top quark decay process, 
both in lowest order and in higher order QCD. As will be seen, this strongly 
reduces the number of independent hclicity amplitudes and makes the results 
quite intuitive. 

2 Helicity vs. Spin Vector Description 

Helicity amplitudes have been considered in many applications-of phenomeno- 
logical importance in high energy physics, like jet productionl!3, nonstandard 
effects in top quark processesli3, and many others. The idea is, first to separate 
a given process into simpler subprocesses and then to explicitly evaluate all 
the possible spin amplitudes for the subprocesses in special Lorentz and Dirac 
frames. The results can afterwards be put together with the help of a master 
formula (to be given below). For example, consider the lowest order hclicity 
amplitudes for top quark decay in the Standard Model. They are given by 

A(h t ,h w ) = -j=u- 1/2 (Pbh^(l ~ 75)u ht (Pt)^hw W 

where g — and h t and hw label the spins for the top quark and the 
W-boson, ht = ±| and hw — 0, ±1. The helicity of the (massless) bottom 
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quark is fixed to be — 1/2 by the V-A nature of the interaction. Thus, in the 
Standard Model there are six amplitudes to be considered, and this number 
remains the same in higher order QCD (after integrating over the gluon degrees 
of freedom). The six amplitudes can be used to define the 36 elements of the 
density matrix 

pih t , h w , h' t , h' w ) := A(h t , h w )A(h' t , h' w )* . (2) 

Note that the amplitudes are only determined up to an overall phase, and 
that this arbitrariness goes away when forming the density matrix elements. 
Unfortunately, higher order QCD corrections cannot be fully calculated on the 
level of amplitudes, but only on the level of the density matrix which contains 
in principle 36 degrees of freedom. As will be discussed below, hermiticity and 
CP invariance reduce this number, but still leave an appreciable set of matrix 
elements. 

The full density matrix is in fact needed in the 'master formula', if one 
considers some combined production and decay process for the top quark. 
Namely, assume that top quarks are produced in some process ab — > tt and 
then decay according to t — > W + b and i — > W~b, where the W's further decay 
to light (massless) fermions, W + — > /1/2 and W~ — > HJa- The cross section 
is then given by the 'master formula' 

cr= I A{h a ,h b ,h u ht)A{h u h w+ )A{h- u h w -)A{h w+ )A(h w -)\ 2 (3) 

EXT INT 

where EXT — h a ,hb denotes the spins of the external particles and INT — 
ht-, hj, h w + , h w - the spins of the internal particles of the process. A(h a , hi,, h tl %) 
are the helicity amplitudes for the production process, A(hi, h w - ) for the de- 
cay of the antitop quark, and A(hw+) and A(hw-) are the amplitudes for the 
decay of the W + and W~ , respectively. Note that before Eq. (||) the W + spin 
has been denoted by hw instead of h]y+ , and for simplicity it will again be 
denoted by hw below. The fact that the formula makes no explicit reference 
to the spins of the massless fermions i=l, 2,3,4, has the same reason as the 
non-appearance of the b-quark spins. Namely, the hf i are fixed by the V-A 
nature of the W decays. Furthermore, if one is not interested in the decay of 
the antitop or of the W's, the corresponding amplitudes and helicities will not 
appear in the above formulas. In that case, the t and/or the W's will be one of 
the external (EXT) particles, whose spins have to be summed over after taking 
the square in Eq. (0). 

Besides neglecting the b-quark mass, it is also a good approximation in 
Eq. (Bh to take the internal particles on-shell, because off-shell contributions 
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are suppressed by powers of the width T t and Tyy . More precisely, one has 

1 n -S(P 2 - M 2 ) + 0( £ ) (4) 



GeV, the width is about 1.5 GeV, so that terms of order -p- can be neglected, in 
particular in higher orders where 0(a s -^r) is of the order of a permille. Similar 



(P 2 — M 2 ) 2 + M 2 T 2 MT ' y M' 

for a particle of mass M and 4-momentum P. For a top quark of mass 175 

r_ 

M 

1 M i 

considerations apply to the W-boson. 

In higher order the narrow width approximation Eq. (|j) is of particular 
use in reducing the number of diagrams to be calculated. The reason is that 
all diagrams where a gluon runs from the production part of the process to the 
decay part, and also the corresponding interference diagrams, give contribu- 
tions which are suppressed by powers of -j^. Therefore, in this approximation 
the process can be really decomposed into a number of building blocks, the 
production block, the t-decay blocks and the W-decay blocks. In the narrow 
width approximation, these blocks are interrelated by spin-indices, but not by 
gluon exchange. 

When one carries out the modulus squared in Eq. (^), it becomes appar- 
ent that in general the full density matrix A(ht, hw)A(h' t , h' w )* is needed to 
calculate the cross section of the decay products. 

I have calculated the six amplitudes A(ht, hw) using the 'chiral represen- 
tation' of 7 matrices, in which 75 =diag(— 1, —1, 1, 1) etc. This representation 
makes calculations with massless fermions (the b-quark in the case at hand) 
quite transparent. Furthermore, I shall present results in two different Lorentz 
frames, in the rest frame of the top quark and the rest frame of the W-boson. 
Both frames have special virtues, so it is worthwhile to study them both. The 
top quark rest frame is of course the natural frame to study top decays, and 
to look at distributions in the energies of the decay products etc. The W rest 
frame has the particular virtue that the amount of longitudinal Ws can be 
read off most easily in this frame. I am quite sure there is a (complicated) 
transformation between the amplitudes in both frames. However, I was not 
able to derive it and, furthermore, found it reasonably convenient to do the 
calculations in both frames separately. 

There is a popular alternative to the helicity formalism, where use is made 
by the fact that the spin of a fermion with 4-momentum P can be described by 
a 'spin vector' S, a pseudo 4-vector which fulfills S 2 — —1 and S-P = 0. In this 
formalism one does not calculate amplitudes but (squared) matrixelements. 
The matrix element for the production of a tt pair with spin vectors St and sj 
has the generic form 

\M\ 2 P ~ tr{j> t + m t )(l + 75/U + m t )(l + 75^) 
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= a + bs t + csi + ds t St 



(5) 



where, for example, d is a tensor with two Lorentz indices and b and d are 
4-vectors. Similarly, the matrix elements for the decay of the top quarks have 
the form 

\M\% = e + fs t \M\% = e + fs- t (6) 
and the full matrix element for production and decay is then given by 

\M\=aee-bfe-cef + dff. (7) 

According to this formula, the cross section will not be just a product of a 
production and of a decay piece, but is given by a suriLof such products. In 
fact, Eq. ^ can be shown to be equivalent to Eq. (^)Eil. 

3 The Method 

The results to be presented were obtained in the helicity formalism. Further- 
more, they concern the oneloop QCD corrections to the lowest order matrix 
elements. The lowest order expressions are usually simple, whereas the oneloop 
expressions, in particular for the case of hard gluons are quite lenghty for ar- 
bitrary spin orientations. Furthermore, the hard gluon contributions cannot 
be treated on the level of amplitudes, because they have to be integrated over 
the gluon's energy and angles. One really has to go to the level of the (spin) 
density matrix, Eq. (Q), to do the phase space integrations. There is, however, 
one circumstance which simplifies the task. This is related to the fact, that the 
oneloop QCD corrections to the total (spin-averaged) width of the top quark 
are known Ei This allows to get rid of the infrared and collinear singularities 
present in the matrix elements, by forming suitable singularity-free combina- 
tions of the spin-dependent and the spin-averaged expressions. The point is 
that the infrared and collinear singularities are 'universal', i.e. independent 
of the spin direction, so that they will drop out in suitable differences. We 
shall discuss our procedure in more detail in the next section, where QCD 
corrections to W — ► qq' are considered as a rather simple warming up exercise. 

4 QCD Corrections to W — ► qq' 

This process is simpler because both outgoing quarks can be considered to be 
massless particles, so that their helicities are fixed by the V — A nature of the 
decay. Therefore, the decay amplitudes depend only on the W-spin and in 
lowest order are given by 

A(0) = sin0 A(±l) = 1±C ° S0 e ±*4> (8 ) 



5 



where 9 and <fi are the (polar and azimuthal) angles between the z-direction 
(defined as the direction of the outgoing quark q) and the direction of the 
W as given in some LAB frame (it may also be considered as the direction 
to which the W-spin points). These formulae hold in the rest frame of the 
W-boson. An overall factor er, X w has been left out in the amplitudes. Note 

2 2 

that the integrated W-width Tw = l""' can be obtained from the trace of 
the density matrix ^2 h A(hw)A(hw)* by multiplying with the square of the 
factor er ^- w and by dividing by the well-known 16-71771^0. 

Using Eq. (g|), the corresponding density matrix Dhw,h' w = A(hw)A(h' w )* 
can easily be calculated in lowest order to yield 

£>'°(0,0) = sin 2 6> 
D lo (±l, ±1) = -(1 + cos 2 6>) ±cosfl 

D lo (+l,-l) = isin 2 6>e 2 ^ 

£>'°(±1,0) = (±cos6>sin6> -sin6> )^ . (9) 

v2 

The underlined terms refer to parity violating effects and cannot be measured 
in hadronic W decays, because the outgoing quark is detected in the form of a 
jet and its flavor and charge cannot be identified. The functions D(hw,h' w ) 
are sometimes called the decay functions of the W. An upper index lo has 
been introduced in Eq. ([)]) in order to make clear that these are the Born level 
contributions to the density matrix. The aim is then to calculate higher order 
corrections to the decay functions/density matrix, in the form 

D(h w , h' w ) = D lo (h w ,h' w ) + ^D ho (h w , h' w ) (10) 

To accomplish this calculation, use was made of the well-known higher or- 
der result for the spin averaged decay function, the 'trace of the spin density 
matrix', which in our normalization is given by 

Dtotal = D\° otal + ^Dt° tal = 2(1 + 2l) . (11) 

The point to notice is that one has 

D(h w ,h> w ) _D l °(h w ,h' w ) a s D h °{h Wl h' w )D\l tal ~ D l °{h w ,h' w )D l Z al 



Dtotal D\° Qtal tt [D\° otal ) 



(12) 
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where the difference D ho (h w , h' w )D l ° otal - D lo (h w ,h' w )D^ tal is completely 
free of singularities separately for hard and virtual gluons. In fact, the diagrams 
with virtual gluon exchange contribute nothing to the ratio D ^^' h w) 
decay, where we shall proceed similarly, the virtual diagrams will contribute, 
but only a finite amount.) The real gluon processes W — > qq'g can be explicitly 
seen to give a finite contribution to the above difference, i.e. the result is finite 
for Eg — > and 9 g —* 0, where-Eg and 9 g are the gluon energy and angle with 
respect to the quark directionlij. The integration over E g and 8 g is therefore 
straightforward and-one obtains the corrections to the W decay functions in a 
very compact form til 

D(h w , h' w ) = (1 + - 0.975^ )[D l °(h w , h' w ) + 0.653^ w ,y 1 . (13) 

TT TT TT w 

This representation can only be obtained after neglecting the irrelevant parity 
violating terms. Furthermore, in higher orders the angles 9 and <fi have been 
defined as to refer to the thrust- instead of the quark-momentum dire cti on. 

These results have been applied to W-pair production at LEP2LL3 and 
NLCLL3. According to a master formula similar to (|3|), the cross section for 
W-pair production and decay in e + e~ annihilation was calculated including 
the higher order corrections Eq. (|l^) and nonstandard contributions. Our 
main motivation in studying this cross section was twofold: 

• first of all we wanted to know how QCD corrections to angular correla- 
tions of W decay products differ from the naive expection of a constant 
K-factor ~ f + Ss.. We found that depending on the kinematic point, 
the deviations from a constant K-factor can be appreciable (of the or- 
der of a few percent). Unfortunately, at LEP2 with its few thousand W 
events these effects are just at the edge to become visible. The situation 
is different at NLC with its larger statistics, where the QCD corrections 
really become relevant. 

• secondly we have proven that QCD correction can mimic the presence 
of nonstandard physics. As has been shown by Monte Carlo studies Ej, 
NLC is sensitive to nonstandard couplings as small as 1CU 3 . This is well 
below the magnitude of the QCD effects induced by Eq. ( |l3| ) tB. 

In addition to the QCD corrections presented above, there are off-shell W 
corrections (with and without QCD E!l) which are particularly important at 
LEP2, i.e. near threshold. Unfortunately, there is no space here to discuss 
them in detail. 
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5 Complete Lowest Order Analysis of the Spin Density Matrix for 

t -► bW 

The decay of the top quark has some extremely interesting physics features, 
which are related to the spin decomposition of the matrix element. In partic- 
ular, it is well known that, due to the large top quark mass value, t-decay is 
dominated by longitudinal W's. In fact, the total width of the W is given by 

r t = r L + r T = - ^|) 2 (i + 2^f ) (u) 



8\/27r mf m 



and the ratio of the number of longitudinal over transverse W's is given by 

2 

Tl/IY = 2 ™?i . Note that the QCD corrections to T t , Tl and Tt are known 

00 to be about -9% resp. 5%, and that all other correction effects like 
electroweak, b-mass and finite width effects contribute roughly 1%. Note 
further that Tl is the most interesting source of loop corrections to the famous 
Rb value [the partial width T(Z — > bb) measured at LEP1]. The point is 
that the exchange of longitudinal W's between the b-quarks gives rise to the 

2 

celebrated corrections of order 0(Gf mf ) ~ Unfortunately, at LEP1 this 

is only a small loop effect which is of the order of 1%, because it is suppressed 
by a factor j^i- In contrast, in t-decays the longitudinal W's enter as the 
most dominant leading order effect, so that they can be studied much clearer. 

Top quark decay may be looked at in different frames; particularly inter- 
esting are the rest frames of the t-quark or that of the W. In the W rest frame, 
for example, there is a very simple way to experimentally determine the ratio 
Tl/Tt and other spin-dependent observableslii. Of course, the two systems 
are related just by a simple boost. However, the transformation formula be- 
tween the spin-amplitudes in the two systems is quite complicated, so that I 
prefer to give results in the two systems separately. Let's start with the t rest 
system. In order to calculate the amplitudes (|l|), I have choosen the follow- 
ing parametrization of momenta, polarization vectors and spinors. First, the 
momenta: 

p t =(m t ,0) pw = — (/+,0,0,/_) Pb=Pt-pw (15) 

where f± — 1 ± / and / = and the W-boson has been chosen to define 
the z-direction. The most general top quark spinor is given by 

/ \ ; / $ • ^ iS 9.9 a, 

u+i/2(Pt) = vWcos-,sm~e <p , cos -, sin -e v ) 

6 9 

u-i/ 2 {pt) = Vrru{- sin -e~ 4< ^ ', cos -, - sin -e~ 10 , cos -) (16) 
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where 6 and <f> refer to some direction, e.g. to the direction of the top quark 
in some lab-system. One may put = without restriction, because this 
corresponds to defining the y-direction. The possible b- and W-polarizations 
are fixed to be 

u_ 1/2 (p 6 ) = ^(0,0 )V /jL,0) (17) 

and 

e= F i = -^(0,±l,i,0) e = --^(/-,0,0,/ + ). (18) 
This leads to the following amplitudes for t-decays in the top quark rest frame: 



1 1 



M-h 0) = 4 7 Binfe- i * A l °(+-,0) = - 



77 0111 ^ ~ t ^r" J -~7j C0S 2 

4°(-i,+l) = 4°(+i+l) = 

4°(-i,-l) = -^co 8 | 4°(+i-l) = -^sin^ (19) 

where the upper index refers to 'lowest order' and the lower index to top 
quark decay in the t rest frame. A universal spin independent coefficient 
c o = \/~f- has been left out in all the amplitudes. Note that the amplitudes 
are only determined up to an overall phase, and that this arbitrariness goes 
away when forming the density matrix elements. Since the amplitudes are 
explicitly given, it is straightforward to obtain the density matrix in lowest 
order. Its trace is easily obtained from the above expressions to be 

Y, 4°(ht, h w )A l °(h t , h w )* = cl 1 -^- (20) 

and one can reproduce from this the total width of the top quark, Eq. ( |l4|) by 
dividing by the phase space factor 16 {. mt El- 
There are lots of other combinations of density matrix elements to describe 
interesting physics. For example, the above mentioned ratio jrp is obtained as 

r, = Z ht = ±hhw=0 \A l t °(h u h w )\* = j_ 

Let us now repeat the same analysis in the W rest frame. This time I 
chose to define the z-direction by the direction of the top quark, i.e. 

Pw = {m w ,0) p t = ^-^(/+,0,0,/_) Pb = Pt-pw- (22) 
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The top and bottom quark spinors are then fixed as 



u +1/2 (pt) = V^(/ 1/4 , 0, /" 1/4 , 0) u_ 1/a (pt) = VrnT(0, Z" 1 / 4 , 0, / 1/4 ) 

(23) 

and 

«-i/a(P6) = VnH(0,0,0, J (24) 



whereas the W polarization direction is arbitrary: 

e_i = t=(0, cos </>cos — z s'mcf), sin^cos + i cos 0, — sin 0) 

e+i = —e*_ x 

e = — (0, sin cos </>, sin sin cos 0) . (25) 

The angles and ^ refer to some arbitrary direction, e.g. to the direction of 
the W-boson in some lab-system. Note that although I am using the same 
symbols and as in Eqs. ( |l9| ) and (@), the meaning of these angles is 
completely different in the three cases. As before, one may in principle put 
= without restriction, because this corresponds to defining the y-direction. 
One is lead to the following amplitudes for t-decays in the W rest frame: 

^w{~,0) = ^f cosd A w(+^ 0) = -sm0e* 
A&(-i +1) = -Lain for* Afc(+± +1) = -L(l + cos0) 
<(-i -1) = --1= sin^ -1) = - cos0)e^ (26) 



Again, the universal coefficient Co = fj^- y /- has been left out in all the ampli- 
tudes. From the trace of the corresponding density matrix one can reconstruct 
the total width of the top quark, just as in (pfj|). However, the ratio ( |20| ) can 
only be obtained for 9 — 0, because otherwise the notion of 'longitudinal' does 
not refer to the heavy quark direction. 



6 QCD Corrections to the Spin Density Matrix of t — ► bW from the 
Exchange of Virtual Gluons 

The oneloop QCD Corrections to t-decay are somewhat more complicated 
than to W-decay because even neglecting the b-mass there is one more mass 
parameter involved. This is despite the fact that the Feynman diagrams needed 
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are exactly the same as for W-decay (with the directions of the W and one 
of the quarks interchanged). Namely, there are the 'virtual gluon' vertex and 
self-energy diagrams and two 'real' diagrams with gluon emission from one of 
the quark legs. As discussed in the introduction, there are 36 density matrix 
elements (||) which will be considered in the normalized form 

(h h h> h' \- P{ h t,h w ,K,ti w ) 

Ptotal 

because this helps to cancel universal contributions, in analogy to the case of 
W decay, Eq. @. Ptotai = J^ht h w \p(ht, hw, fh, hw)\ is defined to be the 
trace of the density matrix and related to the total width Tt of the top quark 
as discussed above. 

Let us start the discussion with the virtual contributions, because they can 
easily be.-abtained from the corrections to the V-A vertex calculated a long 
time agoEy 

r„(t - bW) = --±-{H 7 ^(l 75 ) + ad H + ™^\{l + 75 )} (28) 

where the known function H = 1 + 0(a s ) 'renormalizes' the V-A structure, 
and contains all the ultraviolet, infrared and collinear singularities, and 

H + = (29) 

is a regular function of f — rrv^/rnj. Note that the appropriate expansion 
parameter in the case of top quark decay is ay = — Cf^- It should further 
be noted that the contribution of H to the normalized density matrix p n0 rm 
vanishes because it cancels between the numerator and denominator in Eq. 
( ^7| ) . The argument works in the same way as was discussed in the case of the 
W-boson. The only difference is that now a finite contribution ~ H + from the 
G^v term survives due the nonvanishing (top-)quark mass. In turn, one may 
conclude that the contribution from real gluon emission to the normalized 
density matrix is finite, too, because of the Lee Nauenberg theorem, which 
says that any such singularity cancels between real and virtual corrections. 
One may write down the contribution ~ H + to the amplitudes in the form 

A t (h t ,h w ) = A\°{h u h w ) + a d A+ {h u h w )H+ (30) 



where the lowest order amplitudes A ° were given in (|19|) and the higher order 
amplitudes Af originating from the a^ u term are given by 

At(-hO) = -#sinfe-^ A+(+i,0) = ^ cos °- 
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A+(-|,+l) = Q A+(+i+l) = 

A+ {-I -I) = ^ cos § -1) = sin ^ (31) 



Again, the universal coefficient Co = ff^v/- nas been left out in all the 
amplitudes. These results imply a higher order contribution due to the 
term on the level of the density matrix p{h tl hyy, h' tl h' w ). Including the lowest 
order piece it will be of the form 

A[°{h u h w )A l t °(h' t , h' w )* + a d H+ [A l t °(h t) h w )A+(ti t , h' w )* 

+A+(h u h w )A l °(h' t ,h' w y] (32) 

Note that to calculate p n0 rm, one has to divide afterwards by ptotai, whose 
contribution from the term (including lo) can be calculated to be 

£ \A l °(h u h w )\ 2 + ^-H + ]T [Al°(h t ,h w )A+(h t ,h w )* 

ht,hw ht,hw 

+A+(h u h w )A l °(h t ,h w y] = cgjii^- - 3a d H + } (33) 

This result was obtained by summing over all spin configurations and using 
the explicit representation of the amplitudes given above. It shows explicitly 
that the contribution from virtual gluon exchange to p norm is completely under 
control. 

The amplitudes corresponding to the a^ v term in the vertex may also be 
evaluated in the W-rest frame 

A+(-±,0) = -^cosfl O) = ~sin0e^ (34) 

A+ (-±,+1) = -^ainfcr* A+ (+±, +1) = --L(l + C os0) 
A+(-§,-l) = ^Bin0e* AU+\, -l) = -^(l-cos0)e 2 ^ 



and analgous relations as (|32|) and (|33|) apply. As discussed before, the angles 
9 and </> have completely different meanings in the t and W rest frames. 



7 QCD Corrections to the Spin Density Matrix of t — > from 
Real Gluon Emission 

The amplitudes and density matrix for real gluon emission are the most difficult 
part of the higher order calculation and are also the most difficult to document. 
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The point is not only that the amplitudes are quite complicated expressions, 
but also that one has to deal with the phase space integration over the real 
gluon's degrees of freedom on the level of the spin density matrix. Furthermore, 
the cancellations of the singularities in p norm in Eq. (|27]) require a subtle 
understanding of the interplay between lowest order and first order QCD. All 
this enforces the use of an algebraic computer program like FORM, REDUCE 
or MATHEMATICA to handle the long and complicated expressions. After 
calculating the spin density matrix, the real gluon's degrees of freedom have 
to be integrated over. This integration can be done without regularization, 
because according to the last section the correction to p nor m from real gluons 
is finite and of the generic form, cf. Eq. @, 

ho lo _ lo ho 
(h h h' h' \ J° i ^ ' t 1 total r rtotal (1K\ 

Pnorm,{nt, tlW , "f, tl w ) = P norm + O-d fTo \2 ( 6b > 

\P total) 

where ay — —Cf^- The dependence on the gluon spin is not considered, 
because it cannot be determined experimentally. Accordingly, the gluon po- 
larizations have been summed in the standard fashion. 

To be more explicit, let us write the 4-momenta relevant for the process 
t — > Wbg in the rest system of the top quark as 

p t = (m t ,0) pw = -jr- f+xw (1,0,0, fiw) Pb=Pt-pw~P g 



-^-f—X g (l, sin 9 g cos (f> g , sm9 g sin <j> g , cos8 g ) (36) 



where 



and where 9 g is given in terms of the other variables according to 

f-Xg - /+(1 - x w ) = -f + X W f-Xg(l - (3w cos 8 g ) (38) 

Note that the meaning of /, /_, /+ etc is as in lowest order. The integrations 
over the gluon's degrees of freedom then have to be done with the phase space 

dPS 3 (t -> Wbg) = ^ff + f- f dx 9 d(l - XW ) %L (39) 

4 Jo Jo Jo 27T 

The top quark spinors may be taken as in Eq. (|l^), with <j) — because 
the y-z plane has not yet been specified. However, the W polarization vectors 
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are different from the lowest order expressions, Eq. (|18|), because the W- 
momentum has changed. More precisely, the transverse polarization vectors 
are left unchanged, but the longitudinal polarization vector now reads 

£o = --^jif+xwPw, 0,0, f+x w ) (40) 

The h]y = —1/2 spinor for the b-quark changes, too, but for the density matrix 
one needs only the combination u_i/2W-i/2 = ^(1 ~ 75 )P- 

I have calculated all the 6 x 6 = 36 spin density matrix elements using 
these parametrizations 



(1- 




I -- 

z 


i0) = c 


_ (1 + 0.188 a,;) - 0.0246 c+ a d 


(41) 


(1- 


l~ ^f^Pnorrn (~F 


1 

z 


Is, 

-, 0) = c+ (1 + 0.188 ay) - 0.0246 c_ a d 


(42) 


(1- 


l~ ^f^Pnorrni^ 


1 

0,0,+ 

2 




-s (1 + 0.212 ay) 


(43) 


(1- 




1 

-2'°'- 


i,o,- 


■s (1 + 0.212 ay) 


(44) 


(1- 


l~ ^f^Pnormi^ - 


1 

-2'" 1 ' 




= 2/c+ (1 - 0.236 ay) - 0.00751 c_ 


a^45) 


(1- 


l~ ^'f^)Pnor7n( K ~^ 


1 




= 2/c_ (1 - 0.236 ay) - 0.00751 c+ a£46) 


(1- 


l~ ^f^Pnormi^ - 


1 


^.-') 


= 2/s (l-0.216a d ) 


(47) 


(1- 


h 2/)p norm (4 


1 

2' ' 


-i,-D 


= 2/ « (1-0.216 a*) 


(48) 


(1- 








= -0.00587 c+ a d - 0.0518 c_ a d 


(49) 


(1- 


^ ^f^Pnorm (~F 






= -0.00587 c_ a d - 0.0518 c+ a d 


(50) 


(1- 


^ ^f^Pnorm (~ 






= 0.0460 s a d 


(51) 


(1- 


h 2/)p norm (4 




-> 


= 0.0460 s «d 


(52) 


(1- 


^ ^f^Pnorm ( — 


1 

-2'°'- 




-0.0150 s a d 


(53) 


(1- 


l~ 2,/*)/?norm ( — 


1 

-2" +1 " 


-i,o,- 


-0.0150 s a d 


(54) 


(1- 


^ ^f^Pnorm ("1" 


|o,+ 




+0.0150 s a d 


(55) 



14 



1 1 

l + 2/)p„(+-,+l,+-,0) = 


+0.0150 s a d 


(56) 


l + 2/K orm (--,0,+-,+l) = 


0.0150 c+ a d 


(57) 


1 1 

l + 2/)p„(+-,+l,--,0) = 


0.0150 c+ a d 


(58) 


l + 2/)p norm (+-,0,--,+l) = 


-0.0150 c_ a d 


(59) 


1 1 

l + 2/K orm (--,+l,+-,0) = 


-0.0150 c_ a d 


(60) 


l + 2/)p„(--,0,--,-l) = 


-so y/2f(l + 0.00600 a d ) 


(61) 


1 1 
1 1 

1 + 2/)p norm (--, -1, -2'°) = 


-50^(1 + 0.00600 ay) 


(62) 


l + 2/) / w m (+-,0,+-,-l) = 


+50^27(1 + 0.00600^) 


(63) 


1 1 

l + 2j>„(+-,-l,+-,0) = 


+s oy /2f (1 + 0.00600 a d ) 


(64) 


l + 2/)p„(--,0,+-,-l) = 


-c_ V27(l + 0.00600 a d ) 


(65) 


1 1 

l + 2j>„(+-,-l,--,0) = 


-c_ ^2/(1 + 0.00600 ad) 


(66) 


l + 2/)p„(+-,0,--,-l) = 


+c+ v/2/ (1 + 0.00600 a d ) 


(67) 


1 1 

l + 2j>„(--,-l,+-,0) = 


+c+^2f (1 + 0.00600 a d ) 


(68) 


1 + 2f)p norm {~, + 1, -i, -1) 


= 


(69) 


l + 2/)p„(-i -1,-1 +1) 


= 


(70) 


1 + 2/)/w m (+i +1,+^, -1) 


= 


(71) 


l + 2/) /9 „ orm (+i +1) 


= 


(72) 


l + 2/)p„(+i +1,-1 -1) 


= 


(73) 


l + 2/)p„(-^ -1,+1 +1) 


= 


(74) 


l + 2/) / cw m (- -,+!,+-,-!) 


= 


(75) 
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(1 + 2/) /Wm (+i -1, -I +1) = (76) 

where c± = ^(1 ± cos 9) and sq = | sin 6* and where the numerical coefficients 
in order ay = — Cf§J have been obtained with nit — 175 GeV, for which 
/ = 0.21. With the help of my intergation programs I have shown that the 
m t dependence of these coefficients is in all cases moderate. For example, 
the coefficient 0.188 in Eq. pi] ) depends on = -y= in the way depicted 

in Figure 1. The m* dependence of the other independent QCD coefficients 
(denoted by -0.0246, 0.0212, -0.0150, 0.00600, -0.0518, -0.00587, 0.0460, 
—0.216, —0.236 and —0.00751, i.e. by their values at m t — 175 GeV) are given 
in the figures that follow. Note that only lowest order and real gluons are 
incorporated in Eqs. (^)-(^). The virtual corrections Eqs. (32) and (31) 
from the a^ v term have to be added. The factors 1 + 2/ appearing on the 
left hand side are a relic of the fact that I am presenting the density matrix 
'normalized' to the total width/trace. There are several possibilities to make 
checks on this list. For example, I have checked that 

^ Pnorm(ht,hw,h t ,h\y) = 1 (77) 
ht ,hw 

is true including the oneloop QCD corrections. Furthermore, I have also 
checked that the ratio 

J2h t Pnorm{ht,0,h t ,0) T L 1 ,„ , 

:{l + a s ---) 



Y,h t \Pnorm{ht,+l,h t ,+l) + p norm {h tl -I, h t , T T 2f 



rf 8) 

inM 



reproduces the ratio of longitudinal over transverse W's as calculated int-U in- 
cluding higher order QCD corrections. Finally, there is the check as to the her- 
miticity of the density matrix, p norm (h t , h w ,h' t , h' w ) = pnorm{K-> h 'w h t,h w )*. 
Note that the density matrix is real in the present case, because in the consid- 
ered frame there is not azimuthal dependence. 

I have carried through a second analogous calculation in the rest system 
of the W and obtained the real gluon QCD corrections to the density matrix 
Pnorm{h t , h\y , h' t , h' w ) also in that system. The momenta are now parametrized 
as follows 

pw = (mt,0) p t = -^-jj f + x t (1,0,0, fa) p g =pt-pw-Pb 

TTlf 

Pb = —=f^Xb(l,fimdbCOS(j>b,sm9bSm(j)b,cos9b) (79) 
W J 
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2 2.2 

m(l)/m(W) 



Figure 1: 



-0.015 - 



-0.02 - 



-0.025 - 



-0.03 - 



-0.035 - 




-0.04 1 ' ' ' ' ' ' 1 1 
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Figure 2: 
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Figure 3: 




-0.025 1 ' ' ' ' ' ' 1 1 
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m(t)/m(W) 



Figure 4: 
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-0.05 1 ' ' ' ' ' 1 1 1 

1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 

m(t)/m(W) 



Figure 5: 



-0.035 




-0.085 1 ' ' > > > > 1 1 

1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 

m(t)/m(W) 



Figure 6: 
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Figure 7: 
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Figure 8: 
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Figure 9: 




-0.35 1 1 ' ' ' ' ' 1 1 
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m(t)/m(W) 



Figure 10: 
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2.4 2.6 2.8 3 



where 



Figure 11: 



4/ 

f\4 



and where Ob is given in terms of the other variables according to 
f-X b + /+(! - x t ) = ^jf + x t f-Xb(l - (it cosd b ) . 



(80) 



(81) 



The integrations over the gluon's degrees of freedom are encoded as integrations 
over xt, Xb and <pb and have to be done with the phase space 



dPS 3 (t^Wbg) 



d(pb 
~27 ' 

(82) 

About the effects of gluon emission on polarization: Since the W— momentum 
is unchanged as compared to lowest order, the form of the W polarization 
vectors remains as in (|25|). However, the top and bottom spinors are modified. 
They now read 



"+l/2(Pt) = % /mt(a_,0,a+,0) 



-i/a(pt) = y/rrH(0,a-,Q,a+) (83) 
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with 

a± = ] j^VT±Ft (84) 

and 

«-i/2(p&) = v^H(0.0,&!,&2) (85) 

where b\ and b\ are given indirectly by w_i/2(p&)w_i/ 2 (p{,) = |(1 — 75)^- The 
36 elements of the normalized density matrix obtained in this frame are given 

by 

(1 + 2f) Pnorm {~,0,~, 0) = cos 2 9 + (-0.248 z_ + 0.188 z+) a d (86) 
(1 + 2/)p„(+l 0, +1 0) = / sin 2 + (-0.080 z_ - 0.0246 z+) a d (87) 
(1 + 2/)p„(-l 0, +1 0) = - v7 sintf cos (1 - 0.0171 a d ) 

(l + 2/)p„(+l 0,-^,0) = -^/J sin cos 0(1- 0.0171 a d ) (89) 

(1 + 2/)p„(-l -1, -i -1) = l - sin 2 (90) 

+ (0.188 z_ - 0.0296 z+ - 0.0222 cos 6) a d (91) 

(1 + 2f)p norm {+ 1 -, -1, +i -1) = |(1 - costf) 2 (92) 

+ (-0.0246 z_ - 0.0522 z+ + 0.0464 cos 9) a d (93) 

1 

= ^ [ - sin (1 + 0.0640 ay) + sin (9 cos (9 (1 - 0.0374 ay)] (95) 

(1 + 2/)p„(-i +1, -i +1) = i sin 2 (96) 

+ (0.188 z_ - 0.0296 z+ + 0.0222 cos 6) a d (97) 

(1 + 2f)p norm {+ 1 -, +1, +i +1) = L{\ + cos6) 2 (98) 

+ (-0.0246 z_ - 0.0522 z+ - 0.0464 cos<9) a d (99) 

(1 + 2.f) /Wm (+i +1, ~, +1) = (1 + 2/)p„(-i +1, +1 +1) (100) 

= ^ [ sin 6» (1 + 0.0293 a d ) + sin (9 cos 9 (1 - 0.0171 a d )\ (101) 

(1 + 2/)p„(-i 0, ~, +1) = (1 + 2/)p„(-i +1, -1 0) (102) 
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(1 + 2/)p„(+l -1, -1 -1) = (1 + 2f) Pnorm (~, -1, +1 -1) (94) 



= sin 9 cos 9(1 + 0.218 a d ) - 0.0157 a d sin (103) 
v2 

(1 + 2/)p„(+i 0, +i +1) = (1 + 2/)p„(+i +1, +1 0) (104) 

= — y= [sin6»(l - 0.222 ay) + sin 6 cos 6» (1 - 0.1325a d )] (105) 
v2 

(1 + 2/)p„(-i 0, +1 +1) = (1 + 2/)p„(+i +1, -1 0) (106) 

| cos 0(1 + cos 6») + (0.0130 z_ + 0.00196 z+ - 0.00196 cos 0) ay(107) 

(1 + 2/)p„(+l 0, -1 +1) = (1 + 2/)p„(-i +1, +1 0) (108) 

| sin 2 (9 + (0.00358 z_ - 0.00750 z+ + 0.00750 cos 0) a d (109) 

(1 + 2/)p„(-l 0, -1 -1) = (1 + 2/)p„(-l -1,-^,0) (110) 

= — ^ sine cos 0(1 + 0.218 a d ) -0.0157 a d sin0 (111) 
v2 

(1 + 2/)p„(+l 0, +1 -1) = (1 + 2/)p„(+i -1, +1 0) (112) 

= — y= [sin 0(1 - 0.222 ay) - sin 6 cos (1 - 0.1325 ay)] (113) 
v2 

(1 + 2/)p„(-l 0, +1 -1) = (1 + 2/)p„(+l -1, -1 0) (114) 

^ cos 0(1- cos 6») + (-0.0130 z_ - 0.00196 z+ + 0.00196 cos 9) 4J15) 

(1 + 2/)p„(+i 0, -1 -1) = (1 + 2/)p„(-i -1, +i 0) (116) 

sin 2 + (-0.00358 z_ + 0.00750 z+ + 0.00750 cos 9) a d (117) 

(1 + 2/)p„(-i +1, -i -1) = (1 + 2/)p„(-i -1, -1 +1) (118) 

= -z_ (1 + 0.218 ay) (119) 

(1 + 2/)p„(+i +1, +i -1) = (1 + 2/)p„(+i -1, +1 +1) (120) 

= /*_ (1 -0.1325 ad) (121) 

(1 + 2/)p„(-i +1, +i -1) = (1 + 2/)p„(+l -1, -1 +1) (122) 
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2 
1 



(l + 2/)/w m (+-,+l, 



-1) 



'(1 



ysd) (1 - 0.0171 a d 
= (1 + 2/)p norm (- 
cos 6») (1 - 0.0171 ay) 



1 

2'' 



1. 



1 



-1) 



(123) 
(124) 
(125) 



sin 6>(1 
1 

2' ' 2'" 
2 

where z± = |(1 ± cos 2 ^). The QCD coefficients have again been obtained 
by numerical integration with rat = 175 GeV. The m* dependence of these 
coefficients is again moderate. In fact, it can be shown that the coefficients 
in the rest frame of the top quark, Eqs. (|4l|)-([76|), and of the W boson, 
Eqs. (p^)— (125), are related. I was not able to derive a general formula, but 
I have found, for example, that -0.0518[Eq.@]= -0.0296[Eq. (j9l|)1 -0.0222 
[Eq.©] and -0.00587[Eq.©]= -0.0522[Eq^|)]+0.0464 [Eq.@]. Corre- 
sponding equalities are true for all other values of mt, i.e. they hold for the 
coefficients in general. There are some other relations which I do not want to 
quote here. 



8 Summary 

In this report I have summarized a recent new calculation of a complete spin 
analysis of the Standard Model top quark decay including higher order QCD 
corrections. The QCD corrections to the 'normalized' density matrix are in 
general quite small, of the order of 1%, in particular for the real gluon con- 
tribution in the rest frame of the top quark, cf. Eqs. (pl])-([76|) and of the 
W, cf. Eqs. ( |86|) ( |125| ). The contribution from virtual gluons is somewhat 
larger, cf. Eq. (|29|). Note that the relative magnitude of the QCD corrections 
can be very large - in all cases, where the lowest order contribution vanishes, 
like Eqs. (pi|K(|60|). In other cases, the symmetry requirement of CP gives 
vanishing matrix elements beyond the leading order, cf. Eqs. (p9|)-([76[). 

Complete results including azimuthal dependence, numerical analysis and 
physical applications have not been included here. However, I plan to write a 
long article with J. Korner and his groupL 2 ], in which not only this, but also 
analytical formulae for all the QCD coefficients will be given. 
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